Quantum storage based on the control field angular scanning 
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Continuous change of the propagation direction of a classical control field in the process of its 
off-resonant Raman interaction with a weak signal field in a three-level atomic medium is suggested 
for quantum storage of a single-photon wave packet. It is shown that due to phase matching 
condition such an angular control allows one to reversibly map the single-photon wave packet to 
the Raman spatial coherence grating. Thus, quantum storage and retrieval can be realized without 
using inhomogeneous broadening of the atomic transitions or manipulating the amplitude of the 
control field. Under some conditions the proposed scheme proves to be mathematically analogous to 
the quantum storage scheme based on controlled reversible inhomogeneous broadening of the atomic 
states. 
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I. INTRODUCTION 

Manipulating single-photon states of the electromag- 
netic field is an important part of implementation of 
various quantum information protocols. In particular, 
storage and retrieval of single-photon wave packets is 
of great importance. Optical quantum memory [TH3] 
is at the heart of linear optical quantum computing [4 
and long-distance quantum cryptography with quantum 
repeaters fS. It also may be useful for making her- 
alded single-photon sources deterministic [6 j. An effi- 
cient storage and retrieval of single photons has been 
demonstrated recently in gases [7] and rare-earth-ion 
doped solids [8J using controlled reversible inhomoge- 
neous broadening (CRIB) of resonant atomic transitions. 
Significant progress has also been achieved in demonstra- 
tion of quantum storage using atomic frequency comb 
[9HI4] , electromagnetically induced transparency [T5HT8] 
and off-resonance Raman interaction p~9j [20] . 

It is usually assumed that to store and recall optical 
pulses one needs an inhomogeneous broadened atomic 
transition (tailored or controlled) or a modulated control 
field amplitude matched an input pulse (see reviews pj[3] 
and recent experiments mentioned above). In the present 
work, we develop another approach which requires nei- 
ther inhomogeneous broadening nor temporal modula- 
tion of the control field amplitude, but resorts to con- 
tinuous phase-matching control in an extended resonant 
medium. We consider off-resonant Raman interaction of 
a single-photon wave packet and a classical control field in 
a three-level atomic medium. Under such conditions the 
phase-matching control can be achieved by modulating 
refractive index of the resonant medium [21] or by mod- 
ulating the direction of propagation of the control field, 
which is studied in the present work. Another possibility 



is to use a frequency chirp of the control field, which re- 
quires synchronous modulation of the atomic transition 
frequency in order to keep zero two-photon detuning [22] . 
In any case, a continuous change of the wave vector of the 
control field during the interaction leads to the mapping 
of a single photon state to a superposition of atomic col- 
lective excitations (spin waves) with different wave vec- 
tors and vice versa. In comparison with [21], where the 
cavity-model of quantum memory was considered, here 
we discuss a free-space model. We show that under some 
conditions the proposed scheme proves to be mathemati- 
cally analogous to the longitudinal CRIB-based quantum 
storage, just as a two- level scheme with refractive index 
control does [23 . The angular control scheme, contrary 
to CRIB, allows one to use materials which cannot be 
controlled by external dc electric or magnetic fields. Gen- 
erally speaking, the angular control scheme is also much 
easier for implementation than the refractive index mod- 
ulation. The problem for the latter is that it needs to 
be achieved without modulation of atomic levels. It may 
be done in some specific materials like Tm:LiNb03 [23] - 
But in general case such scenario is rather difficult to 
implement as it was discussed in [24] . 

The paper is organized as follows. In Sec. [TT] the model 
of the quantum memory is presented. In Sec. 
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analyze the storage and retrieval of single-photon wave 
packets with transverse propagation of the control field. 
In Sec. |IV[ a more general situation is considered in the 
limit of a large Fresnel number of interaction volume. In 
Sec. [V| we consider some implementation issues. Sec. [VI 
concludes the paper with final remarks. 



II. THE MODEL AND BASIC EQUATIONS 

We consider a system of identical three-level atoms in- 
teracting with a weak quantum field (single-photon wave 
packet) to be stored and with a strong classical control 
field. The atoms have a A-type level structure, and the 
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fields are Raman resonant to the lowest (spin) transition, 
see Fig. []Ja). We assume that the atoms are stationary- 
like impurities embedded into a solid-state material or 
cold atoms in an optical lattice. The sample is approx- 
imated by a parallelepiped with cross section L x x L y 
and length L. The coordinate system is originated at the 
center of the medium. The single-photon wave packet 
is x polarized and propagates in z direction. It is as- 
sumed that the transverse spatial profile of this quantum 
field completely goes into the sample cross section. The 
control field is supposed to be polarized along axis y and 
propagated in any direction in (x, z)-plane. The duration 
of storage and retrieval processes is equal to T. During 
the storage, t G (— T, 0), the propagation direction of 
the control field changes so that at different moments, 
the signal field interacts with the control field of differ- 
ent directions, thus creates coherence between |1) and 
13) (spin wave) with different spin wave vectors, see Fig. 
[lib). During retrieval, t G (0, T), the interaction between 
the spin wave and the control field reconstructs the sig- 
nal field, as illustrated in Fig. [TJc). The storage time 
between the end of the storage and the beginning of the 
retrieval is subjected a free decay of the spin wave, which 
will be neglected through out this paper. 




FIG. 1. (Color online) (a) Energy diagram of the Raman 
interaction in a three-level A system; (b) During storage, the 
temporal profile of the signal field is mapped into the spin 
wave distributed over different spin wave vectors; (c) During 
forward retrieval, the spin wave profile is mapped back into 
the output signal field. 

The signal field corresponding to the single-photon 
wave packet of average frequency uj s and wave vector k s 
is considered in the paraxial approximation and written 



as 

E s (r,t) = ^ x [^a(r,t)e i ^-^+R.c, (1) 
n V 2£ c 

where n and c are the refractive index and phase ve- 
locity of the quantum field inside the medium, a(r, t) is 
the slowly varying annihilation operator. Expression ([I]) 
is written in the approximation of the equal phase and 
group speeds of light in the medium (see, e.g., [25]) which 
is made for simplicity. 

The classical control field is supposed to be a 
monochromatic plane wave subjected to a phase modu- 
lation due to rotation of its wave vector on a small angle. 
Since most of the atoms remain in the ground state dur- 
ing the interaction, we neglect the control field absorption 
so that the field can be written as 

E c (r, t) = E e *(k c r-^+0(r,t)) + c c ^ (2) 

where E is a constant amplitude of the plane wave, k c 
is an average wave vector of k c (t) during rotation, and 
0(r, t) is a phase shift due to the rotation, which is spec- 
ified later in this section. 

The collective atomic operators are defined as the 
mean values of the single-atom operators 

cr mn (r,t) = i ^K-)(n,-|^ 3 )(r - r,) , (3) 
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where N is the atomic number density, which is supposed 
to be constant in space, and \rij) is the nth state (n = 
1, 2, 3) of jth atom with the energy hu n (lji = < us < 
CJ2). The slowly varying amplitude describing coherence 
on the Raman transition, s(r, t), is then introduced as 

ai 3 (r, t) = s(r, t) e *(k s -k c )r-;(^-^ c )t (4) 

Off-resonant Raman interaction is described by the fol- 
lowing equations: 

f d 1 d \ . A i . . fN 

+g*Ns{v,t)e i ^ r ' t \ (5) 
|U(r, t) = (-7 + iS) s(r, t) - ga(r, t) , (6) 

where A ± = £ + g = \ >h£ is the coupling 

constant between the atoms and the weak quantized field, 
Q = ds2Eo/his the Rabi frequency of the classical control 
field, dij is the dipole moment of the transition between 
\i) and A = U2 — uj s is one-photon detuning, 7 is the 
rate of dephasing of the spin coherence, which in the gen- 
eral case includes both homogeneous and inhomogeneous 
broadening of the Raman transition, and 5 is two-photon 
detuning. When writing Eqs. ([5| and |6|, it is assumed 
that 1) the time of propagation of photons through the 
system L/c is negligibly short compared to the evolu- 
tion time of the slowly time-varying envelopes; 2) all the 
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atoms are initially in the ground state |1) and most of 
them remain in the ground state so that Langevin noise 
atomic operators are not included; 3) the frequency shift 
|S7| 2 /A of the Raman transition induced by the coupling 
field is taken into account by redefining field frequency 
uo c ] 4) the refractive index change of the medium due to 
the atoms is incorporated into the value of n. 

It is convenient to expand the slowly varying atomic 
and field operators on the transverse mode basis 



([5| and (|6| in the reciprocal space take the form 



d 19' _ 

dz + ~cdt ]Arnn[z ' t] ~ 



g*NJ2Srn>n>(z,t)B+, n ,^ n (z,t), (15) 



_9 

at 



a ( r > ^ = J2 U mn( x i y)amn(z, t) , (7) 
mn 

y)s mn (z,t), (8) where 



S mn (z, t) = (-7 + iS) S mn (z, t) 
-g A m>n' (Z, t) B-, n , mn (z, t) . 



(16) 



where the transverse mode functions u mn (x, y) satisfy the 
conditions of completeness and orthogonality: 

^2 u *mn( x i V) v>mn(x\ y) = 5(x - x) S(y - y) , (9) 
dxdyu* mn (x,y)u m ' n '{x,y) = 5 mm > 5 nn > . (10) 



In what follows, we use the set of paraxial plane waves 
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(11) 



where p = (x,y), q mn = (2nm/ L x ,2nn/ L y ), m,n e Z, 
and L x1 L y are transverse linear sizes of the sample of a 
rectangular cross section. This definition corresponds to 
periodic boundary conditions in the x, ?/-plane. 

In order to write down the equations in the transverse 
reciprocal space, the function 0(r, t) needs to be specified. 
Let us suppose that the wave vector of the control field is 
rotated around some average direction on a small angle 
so that the induced phase shift may be considered as 
a linear function of the rotation angle. In such linear 
regime, we can write 

0(r, t) = p x (x - x )t + P y (y - y )t + fi z (z - z )t + o (r), 

(12) 

where Pi is the rate of change of the component k c along 
ith axis, and the coordinates #o, 2/o? z o correspond to a 
phase stationary point where the phase of the control 
field remains constant during the rotation. 0o(r) is a 
time-independent phase factor which can be incorporated 
into s(r, t). 

In addition, we extract factors describing the phase due 
to transverse momentum by introducing new variables 



(13) 
(14) 



where w mn (z) = exp 



' 2k, 



<z 



, and z w corresponds 



B t>n>, m n( Z it) = Sinc 



i^Px^ H~ Qm'n' ,x Qmn,x^) -Lx 



smc 



B 



W* m , n , (Z) W mn {z) e ±itf-*-/J-*o-/5.xo-/5 sM) )t _ (17) 



Without rotating the control field the factors 
m'ri mn( z ^) become identity matrices, and each trans- 
verse mode can be considered independently in the frame- 
work of one-dimensional model. However, even in the 
case of rotating control field, the decoupling is also pos- 
sible if (3 x t and f3 y t remain small with respect to the inter- 
val between the modes 2tt/L x and 27i/L yi respectively, 
while p z t being larger than 2tt/L. In such a case, different 
Fourier transverse modes evolve independently on each 
other. The quantum storage can then be considered in 
single-mode approximation, which is the subject of Sec. 
Ill Another simplification is possible in the limit of a 
large Fresnel number F of the excitation volume. In this 
situation, diffraction spreading may be neglected, and 
Eqs. (|5| and (|6| allow one to consider each point (x, y) 
independently. Such a geometrical optics approximation 



is considered in Sec. |TVj Regarding the excitation geom- 
etry, it should be noted that fi x xo -\-(3 y yo-\- fi z zo becomes 
zero when average control wave vector k c is directed from 
stationary phase point (xo,^o 7 ^o) to the center of the 
sample. We will see in Sec. IV that this is the most 
convenient excitation geometry, which needs a zero two- 
photon detuning S and requires no manipulation with S 
during storage and retrieval. 

The main figures of merit those should be calculated 
for the quantum storage are total efficiency 77 and fidelity 
& . Since the input and output faces of the sample are lo- 
cated at the points z = —L/2 and z = L/2, respectively, 
the input and output single-photon wave packets are de- 
scribed by the operators A mn ^ n (t) = A mn (—L/2,t) and 
Amn iOU t(t) = A mn (L/2,t). The total efficiency is defined 
as 



to the position of the beam waist. Then the set of Eqs. 



(18) 
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where 



Nin = Y, dt(Al n . n (t)A 

mn, in (*)>, (19) 

mn 

POO 

N out = / dt(Al^ out (t)A mnjOUt (t^ , (20) 

mn J 

since we assume that the storage process terminates at 
the moment t = 0, while retrieval process begins at this 
moment of time. 

For the time-reversed output pulse the fidelity is usu- 
ally defined as 



where 



X^mn Jo ^ (^mn, out (^)^-mn, in 



out 



(21) 



(22) 



is a measure of pulse preservation independent of 
the total efficiency, and i is the delay that maxi- 
mizes . In what follows, it is the latter quan- 
tity that is used for characterizing quantum stor- 
age together with the efficiency. The definitions of r\ 
and can also be done in a similar way in terms 
of the spatial variables ai n (r, t) and a out (r,t). In 
real space 7V in = / d 2 p j_ OQ dt (a\ n {p,t)a in {p,t)) , and 
^out = Jq 00 dt J d 2 p(al ut (p,t)a out (p,t)}, while = 

Nin N out Io° dtfd 2 p(al n (t-t)a out (p,t)) . 



III. SINGLE-MODE APPROXIMATION 

For simplicity, we assume in what follows that the 
wave vector of the control field is rotated in (x, z)-plane 
around an average polar angle #o, as shown in Fig. |2j 
so that k c = /c c sin#o x + /c c cos #o z ? where k c = 
This also ensures that the control field polarization is 
unchanged during the interaction. In this case, we have 
f3 x = /3cos#o, /3 y = and /3 Z = — f3 sin 6n. The total angle 
of rotation 2A0 during the storage or retrieval process is 
2A0 = f3T/k c = /3TX c /2tt <C 1. 

Now we consider the case when the control field prop- 
agates perpendicular to the signal field, i.e. On = 7r/2, 
which corresponds to a transverse control field. Retrieval 
is done in a forward way. The transverse size of our sam- 
ple is assumed to be smaller than the sample length. In 
this case, as pointed out in Sec. [TTJ it is possible to switch 
between different longitudinal modes without switching 
between transverse modes when rotating the control field. 
The angle of rotation AO should satisfy the following con- 
dition in order to stay within a single transverse mode: 



AO < 



2T A, 
At L x 



where At is the duration of the signal 



pulse. In such a case, Eqs. (15) and (16) are decou- 
pled, and w mn gives only diffraction effect for each of 



the transverse mode of the signal field. Thus it is suffi- 
cient to consider the evolution of a single mode. Let us 



define S„ 



z ° and go to co-moving frame 



r = t — z/c. Since [3 X = and (3 Z 



from Eqs. (15) 



and (16) we have 



r o 



dz A mn (z,T)=g*NS mn (z,T), (23) 
s 'mn( z , T ) = ("7 + 0(z - z j) S mn (z, t) - 



dr 



-gA mn {z,T) . 



(24) 



Then one recognizes that in such a regime, equations 
describing the system are the same as those for longitu- 
dinal CRIB scheme ( which is also referred to as gradient 
echo memory scheme) [26]. In longitudinal CRIB scheme, 
a space dependent absorption line along the medium is 
created. While the signal field propagates through the 
medium, different frequency components get absorbed 
by different absorption lines at different longitudinal po- 
sitions, resulting in a space dependent coherence in the 
sample. During retrieval, such coherence is mapped back 
onto the output signal. We can understand our scheme 
in exactly the same way: Since the rotation of the control 
field yields a factor e l ^ r ^ = e l ^ zt ^ z ~ z °\ thus on the one 
hand, there is a time dependent wave vector fit which 
is responsible for the writing of the spin waves with dif- 
ferent wave vectors and the recording of the temporal 
profile of the signal field. On the other hand, this term 
can be viewed as a space dependent absorption line at 
frequency /3z, which absorbs different frequency compo- 
nents of the signal at different positions of the medium. 
The medium opens an absorption window of width f3L 
(Fig. Ilia)) along longitudinal direction. The absorp- 
tion of the central frequency of the signal field happens 
at z p = zo + S//3. It is necessary for z p to be inside 
the medium, and better at the middle of the sample. 
This means that if zo = 0, i.e. the longitudinal position 
of the phase stationary point of the control field corre- 
sponds to the center of the sample, then the two-photon 
detuning 5 should be equal to zero. Otherwise, one can 
take advantage of the two-photon detuning to shift the 
phase stationary point longitudinally to the medium cen- 
ter. Although the equations of motion can be reduced to 
the same as those of CRIB, the underlying physics of 
our scheme is quite different. Instead of controlling in- 
homogeneous broadening, we realize quantum storage by 
continuously creating spin wave vectors of different val- 
ues to record the temporal information of the incoming 
single-photon wave packet. In other words, we achieve 
the same storage effect as in CRIB via a specific phase 
control of the control field which results in the same ef- 
fect as the spatially depended frequency control of the 
atomic levels. 

In order to have good storage efficiency, the absorption 
window width f3L should cover the input pulse spectrum 
width: f3L > 2^/ At. So 



2A0 > 



At L 



(25) 
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This condition is the same as that resulted from switching 
between different longitudinal modes during the storage 
of the signal pulse. On the other hand, it is known that 
the "optical density" for each spectral component of the 
signal field in Fourier space for longitudinal CRIB [27] is 
2tt \g\ 2 N//3. This quantity needs to be larger than unity, 
so 



2A<9 < T\ ( 



\q\ 2 n. 



(26) 



Thus one needs to keep a balance between the absorp- 
tion window width and the "optical density" . The bigger 
P is, the wider the absorption window opens while the 
smaller the "optical density" becomes, and vice versa. 
Besides, we choose the parameters to avoid significant 
storage phase factor [27] throughout the paper. The re- 
trieval of the signal is done by switching the sign of j3 
(and the sign of J, if a non- vanishing 5 is used to shift 
the absorption center), corresponding to a reversal scan 
of the control field. The retrieval field experiences a phase 
modulation [28], which can strongly decrease the fidelity 



defined as (22), and has to be taken into account be- 



yond conditions (25) and (26). The performance of the 



storage and retrieval is then equivalent to that of longitu- 
dinal CRIB, which has been well studied and understood 
in many papers [28-31. However, it is worth noting that 
a retrieval signal can also be generated without switch- 
ing the scanning direction of the control field (thus the 
sign of /3), as opposed to longitudinal CRIB scheme. Yet 
in this paper, we always perform a reversal scan during 
retrieval. 

We note here that Eqs. ([23]) and ([24]) with z = 



—L/2 are also the equations describing quantum mem- 
ory scheme via refractive index control [21 in free space, 
for which the above discussions are applicable. In this 
case, a non- vanishing two-photon detuning S — PL/ 2 is 
required to shift the position of the absorption line of the 
central frequency of the signal field to the center of the 
medium. 



IV. GEOMETRICAL OPTICS 
APPROXIMATION 

In the previous section, we consider a special case when 
the propagation directions of the signal and the control 
fields are perpendicular to each other. This allows us 
to decouple equations in the reciprocal space and make 
use of single-mode approximation. However, the trans- 
verse excitation requires more power because of the large 
control beam cross section. In this section we consider 
control field propagating at an arbitrary angle with re- 
spect to the z axis. But we make use another approxi- 
mation when the Fresnel number F of interaction volume 
is much bigger than unity. Then the second order trans- 
verse derivative in Eqs. ([5| can be dropped so that each 
point in the transverse plane (x,y) can be considered 
independently. The present case corresponds to geomet- 
rical optics when the signal field is described as a series 



the second order 
2fi-Aj_a(r,t) gives 



of rays propagating along the axis z. The transverse pro- 
file of the signal field is assumed to be smooth enough 
and treated as featurelessness within each of the Fresnel 
zones. 

In the reciprocal space, 
transverse derivative term 
—iq^ nn /(2k s )a rnn (z 1 t). This term describes phase 
shift due to transverse momentum, which can be ignored 
when q^ nn /k s <C n. The latter corresponds to a very 
directional propagation of the signal field when Fresnel 
number of interaction volume is large. However, when 
the control field is rotated and different spin waves are 
created in the medium during the Raman interaction, 
the signal beam can be subject to spreading in transverse 
direction due to the scattering of the control field on 
the spin waves with different wave vectors. Therefore, 
in order to neglect the transverse derivative in the 
equations, we need not only large Fresnel number of 
the interaction volume, but also sufficiently small angle 
of control field rotation AO. The upper limit of AO 
can be estimated in the following way. Consider the 
interaction between signal field of wave vector k s and 
control field of wave vector k c (t) (Fig. [2|. While k c 
changes from k c (ti) (at an angle Oo — AO) to k c (^) 
(at an angle Oq + A0), it continuously creates a set 
of spin wave vectors k s — k c (t) depicted in Fig. [5J 
As a result the signal wave vector can spread over 
an angle cp. For k s ~ fc c , AO <C 1, one can estimate 
(p < 2A0. We need (Aq% + Aq%) L/ (2k s ) <C tt. Since 
Aq x ~ /c s sin(/9 ^ 2A0k Sl Aq y = 0, we arrive eventually 

at the following limit: AO <C A# max = ~ \ VF0 d , 

where 0& is the diffraction angle. It should be noted 
that this inequality can be considered as just a sufficient 
condition for neglecting second derivatives and geomet- 
rical optics approximation, but not necessary for the 
considered quantum storage protocol. 

Now, by defining a new variable s'(r,t) = s(r, t)- 
, e i[i3 x (x-x )t+i3 z (z-z )t] ^ we h ave ^he following equations 

in the co-moving frame: 



r d_ 

dz 



a(r,r) =^*7V 5 / ( r ,r), 



(27) 
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s'(r, r) = {—7 + iS + i[P cosOq (x — xq) 



-PsinOo (z - zo)]K(r,r) - ga(r,r) . (28) 



dr 



Taking the same argument as in Sec. |III[ we arrive at 
the following condition for efficient quantum storage: 



^^<2A<9sin6>o<TA c |#| : 
At L |y| 



N. 



(29) 



Beyond this, we need to confine the value of the average 
polar angle 0q. The central absorption line appears at 
the position z p = „ J n Qq + zq + cot 0q{x — xq). z p should 
be inside the medium, so 



L/2 <z p < L/2. 



(30) 
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F > 1 condition (32) automatically ensures that 9q 3> 



b) 




FIG. 2. (Color online) Diagram of the signal wave vector, 
control wave vector and spin wave vector for (a) forward re- 
trieval and (b) backward retrieval. A writing control wave 
vector changes direction from k c (ti) to k c (to). For t\ — — T, 
#2 = 0, the control wave vector sweeps an angle 2 AO around 
an average polar angle # . When a reading control field is 
applied at time t\ the retrieved signal field spreads over an 
angle ip. In backward retrieval, the spin waves need to be 
flipped to their conjugated components in the time interval 
between the end of the storage and beginning of the retrieval. 



Again, one can engineer the control field to make the lo- 
cation of the phase stationary point (#o, ^o) = (0,0) and 
set two-photon detuning 5 = 0. However, this condition 
can be relaxed in a way that the frequency shift intro- 
duced by zo is canceled by the shift introduced by xo and 
5, namely, zq — xq cot 0q + 5/ (/? sin O ) = 0. So the phase 
stationary point does not need to be exactly at the center 
of the sample, but may be situated anywhere on the line 



xq = zq tan O + 8/ {P cos 6q) 



(31) 



When (5 = 0, this is just the bisector of the rotation. If 
there exists any difficulty to put (xq, zq) on the bisector, 
a non-vanishing 5 can be used to shift the position of 
the absorption lines. In the following we will assume this 
condition is fulfilled. Then (30) gives —L/2 < xcot^o < 
L/2. Replace x by the rac 
R, we have |cot#o| < 
arctan then the condition on the average polar angle 
On can be written as 



ius of the excitation volum 
Define geometry angle g = 



% > g = arctan 



2R 
~L 



(32) 



If (31 ) is not fulfilled, the angle 0q needs to be even bigger. 
Since g ~ F0d and A# max « \\fF0d-, for Fresnel number 



A0, which is the condition we have been using through 
out this paper. 



Forward retrieval 



Under the condition ([31]), the Eqs. ([27]) and ([28]) be- 
come 



( d 

— a(r,r) = g*Ns'(r,T), 
d 

— s'(r,r) = [-7 + i/3(x cos <9 



(33) 



-z sin (9 )]5 / (r, r) - ga(r, r) . (34) 



The term /3x cos 6 ^s' r (r,t) has no positive contribution to 
the storage and retrieval of the signal field. It appears 
as a side effect due to the change of k c on transverse 
direction. For forward retrieval, this term could intro- 
duce some transverse profile to the retrieved signal. This 
transverse profile reveals itself as a transversal space de- 
pendent time shift (ahead or behind) of the forward re- 
trieved signal, as seen in Fig. [3](d). The reason of the 
time delay is as follows: the signal field in the medium 
excites spin wave and evolves in the form of polar iton, 
which is the combination of the signal field and the ex- 
cited spin wave [29]. The polariton, after storage, freezes 
inside the medium at the position of central absorption 
line z p = x cot 0q • It is now clear that the condition ( 30 ) 
is to place the polariton just inside the medium. Other- 
wise if no polariton is created, the signal field can not be 
stored. During forward retrieval, the control field picks 
up the polariton at z p and generates retrieval signal at 
z = L/2. So the group velocity of the signal field reduces 
from c to while z = —L/2 — >> z p , and recovers from 
to c while z = z p — >> L/2. As a result it is obvious that 



if z p is not exactly equal to 0, there must be a time shift 
going along with the forward retrieved signal. This time 
shift is avoided in the case of 0o = tt/2 because there z p 
is independent on x (transversely homogeneous). Even 
if z p is not exactly equal to 0, there would be only a 
time shift as a whole for the recalled signal. However, 
for 0o ^ 7r/2, at x ^ this time shift is unavoidable 
for forward retrieval, and moreover, it must introduce a 
transverse distortion of the signal. For z p not deviating 
much from and 4c \g\ 2 N/(f3 2 L 2 sin 2 O ) > 1, the group 
velocity of the signal field near the ends of the medium 



is [28] v g (z)\ lz±L/ 2\ <L /2 ~ P 2 (z - Zp) 2 / ( 
the delay is due to the propagation between z 
and z = L/2, from dz 
estimated as 



g\ 2 N 



Since 



L/2+2z p 

v g (z)dt, the delay time can be 



t d 



\g\ 2 N 



i 



i 



f3 2 sm 2 \L/2 + z p L/2 - z T 



(35) 



Taking z p = ±RcotOo, one estimates the temporal 
broadening of the forward retrieved signal with respect 



7 





(a) 1- 
0.98/ 
0.96 
0.94 
0.92 
0.9- 



tt/8 



(c) 1- 
0.98' 
0.96 
0.94 
0.92 
0.< 



k/S 



71/4 



3ti/8 



7C/2 



Tt/4 



3ti/8 



7C/2 




0.6 
0.5 



tt/8 



tt/4 3tt/8 



n/2 



FIG. 3. (Color online) Forward retrieval fields for 0o = tt/2 
in (a) & (c) and 0q = ty/9 in (b) & (d). (e) shows the 
signal field to be stored. The vertical axis x/(L x /2) is the 
transverse dimension of the sample normalized by the half 
transverse size L x /2\ the horizontal axis r/T is the local time 
normalized by the duration of the storage (retrieval) process 
T. The color indicates the signal field amplitude, (a, b) 
Retrieval fields for constant AO. The field is distorted for 
small #o because (1): the absolute value of f3 z — — /3sin#o 
becomes smaller such that it approaches the lower limit of 
(29); (2): f3 x — f3 cos 0o becomes large enough to lead a 
spacial-temporal distortion; (c, d) Retrieval fields for con- 
stant AO sin #o • There is transverse distortion of the retrieval 
field due to the non- vanishing /3 X . The figures are gener- 
ated under the following parameters: A s « A c = 1550 nm( 
fiber-optic communication band), \g\ 2 N — 8.3 x 10 10 /(s- m), 
At/T = 1/20, T = 1000 ns, 2R/L X = 1/6, where R is con- 
sidered to be the transverse spatial half width of the input 
signal, L x — 0.6 cm, L — 1 cm, A^ = 8x 10~ 3 rad in (a) & 
(b), = 8 x 10" 3 / sm6, o rad in (c) & (d). 



to the input signal as 



g\ 2 N 



4R cot 6» 



/? 2 sin^o (L/2)-R 2 cot 2 o 



(36) 



which is finite and positive considering the condition (32 ). 



Fig. [4] shows how the efficiency and fidelity depend 
on #o- Although a smaller 0q requires less control field 
power, it gives less fidelity due to the transversal dis- 
tortion for a given AO. Alteratively, it is necessary to 
increase A6 when 0o approaches g from above, in order 
to maintain high fidelity of quantum storage, as shown 
in Fig. [I](c, d). 



B. Backward retrieval 

In the case of forward retrieval, the above transverse 
distortion decreases the fidelity of quantum storage when 



FIG. 4. (Color online) Forward retrieval Efficiency 77 and Fi- 
delity v.s. the average polar angle 0o. (a, b) Efficiency (a) 
and Fidelity (b) for forward retrieval with AO kept as constant; 
(c, d) Efficiency(c) and Fidelity(d) for forward retrieval with 
A#sin#o kept as constant. The parameters are the same as 
described in Fig. [3] 



the average polar angle #0 becomes small. This can be 
avoided in a backward retrieval, wherein the time shift 
at each transverse point is exactly compensated during 
backward propagation. As a result, the fidelity remains 
high for all average angles #0 of interest. However, due to 
the phase matching condition, it is not enough to merely 
switch the propagation direction of the control field to 
the opposite. In general, it requires phase conjugation 
of the spin waves. Under considered geometrical optics 
approximation, it is sufficient to flip the direction of the 
spin wave vector along direction K = k s — k c (£ = 0) 
before the retrieval. This can be done by applying two 
non-colinearly propagating tt pulses one followed by an- 
other [32j|33], both on resonance with level |3) and an- 
other energy level, say, |4). The transition frequency CJ34 
needs to be bigger than e|K|. The directions of the two tt 
pulses are arranged in such a way that the wave vectors 
difference between them is along — K. Another possibil- 
ity is to generate spin wave vectors perpendicular to the 
signal wave vector [34]. In this case, k c and k s needs to 
differ from each other significantly. 

In the case of backward retrieval, the input signal en- 
ters at z = —L/2 and creates polariton at z = z p: while 
the retrieved signal starts from z = z p and emits at 
z = —L/2. As a result, there is no delay of the output 
field for any value of z p G (—L/2, L/2). 

Generally speaking, in either forward or backward re- 
trieval, the scanning of the control field can be imple- 
mented in different ways. The control of the scanning 
order and scanning rates offers the possibility of the ma- 
nipulation of the retrieval pulses, for instance, recalling 
different temporal parts of the input pulses in different 
orders. Such a manipulation in the angular-time do- 
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main may be thought of as counterpart of that in the 
frequency-spatial domain using CRIB [35] . 



V. IMPLEMENTATION ISSUES 

Let us discuss possible experimental conditions under 
which the storage and retrieval of weak optical pulses via 
angular scanning of the control field is possible. 

In Fig. [3] and [4] we demonstrate the performance of 
our quantum memory scheme with signal field frequency 
around fiber-optic communication band. Yet experimen- 
tally it is preferable to store and recall the signal field 
with small total angle of rotation 2A6 and rotation rate 
2A0/T. If we take A s « A c = 700 nm, according to the 
numerical simulations, a single Gaussian pulse of dura- 
tion At = T/20 with T = 500 ns can be stored and 
recalled with high efficiency and fidelity using transverse 
control field if the angle of rotation is of the order of 
10A/L, which gives AO/T ~ 10 3 rad/s for L = 1 cm. 
Such rate of beam deviation can be achieved by com- 
mercial equipment. Regarding promising storage mate- 
rials, two points must be kept in mind in the context 
of the present scheme: 1) off-resonant Raman interac- 
tion needs systems with relatively strong optical transi- 
tions, especially in the case of transverse control field; 
2) exploiting coherent spatial grating needs atoms to be 
stationary in space. From this point of view, a system 
of cold atoms trapped in an optical lattice seems to be 
one of the promising candidates [36j [37] . The influence 
of regular atomic structure on the efficiency of quantum 
memory in such a system was analyzed recently in [38] . 
We note only that the maximum length of the wave vec- 
tor of a spin wave created via non-collinear Raman inter- 
action is limited by the minimum interatomic distance, 
which means that a blue detuned optical lattice should 
be used for the case of transverse control field. Another 
promising material is an ensemble of defect centers in 
diamond such as nitrogen- vacancy centers [39j [40] . The 
existence of lambda type optical transitions in such a sys- 
tem was demonstrated via electromagnetically induced 
transparency [41, and coherent population trapping [42] . 



The effect of inhomogeneous broadening of the Raman 
transition may be removed through the use of a spin-echo 
pulse sequence, and long dephasing times are achieved by 
dynamic decoupling of spin qubits from their local envi- 
ronment [43-45]. Moreover, transfer of quantum states 
between the electron spins and nuclear spins is possible 
[46] , which in combination with dissipative decoupling 
schemes allows long-lived quantum storage [47] . 

VI. CONCLUSION 

We propose a new method to store and retrieve weak 
pulses such as single-photon wave packets based on off- 
resonant Raman interaction. By changing the propaga- 
tion direction of the strong classical control field, the tem- 
poral profile of the signal field is mapped into the spatial 
grating of Raman coherence. If the wave vector of the 
control field is perpendicular to that of the signal field, 
the quantum storage can be described in single-mode ap- 
proximation so that the proposed scheme proves to be 
equivalent to longitudinal CRIB scheme. When the con- 
trol field approaches longitudinal one, in the case of for- 
ward retrieval there maybe additional spatial-temporal 
distortion of the output field with respect to the input 
signal due to the change of the control wave vector on 
transverse direction. Such distortion can be avoided in 
backward retrieval. The proposed scheme has the ad- 
vantage of longitudinal CRIB in that high efficiency can 
be achieved without backward retrieval. Besides unlike 
CRIB, this scheme does not require a direct control of 
atomic levels, thus potentially reduces decoherence in the 
system. The scheme can be implemented in resonant me- 
dia which do not demonstrate linear Stark or Zeeman ef- 
fects and allows one to combine Raman-interaction-based 
and CRIB-based approaches solely in the framework of 
the former. 
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